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This Article presents a three dimensional numerical model investigating thermal performance and hydrodynamics 
features of the confined slot jet impingement using slurry of Nano Encapsulated Phase Change Material (NEPCM) 
as a coolant. The slurry is composed of water as a base fluid and n-octadecane NEPCM particles with mean di- 
ameter of 100nm suspended in it. A single phase fluid approach is employed to model the NEPCM slurry. The 
thermo physical properties of the NEPCM slurry are computed using modern approaches being proposed recently 
and governing equations are solved with a commercial Finite Volume based code. The effects of jet Reynolds 
number varying from 100 to 600 and particle volume fraction ranging from 0% to 28% are considered. The com- 
puted results are validated by comparing Nusselt number values at stagnation point with the previously published 
results with water as working fluid. It was found that adding NEPCM to the base fluid results with considerable 
amount of heat transfer enhancement.The highest values of heat transfer coefficients are observed at H/W=4 and 
C,,=0.28. However, due to the higher viscosity of slurry compared with the base fluid, the slurry can produce 


drastic increase in pressure drop of the system that increases with NEPCM particle loading and jet Reynolds 


number. 
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Introduction 


Heat transfer enhancement in jet impingement process 
is one of the most reliable technique for high heat flux 
removal from heated surfaces in many engineering and 
industrial applications such as cooling of gas turbine 
blades and electronic chips, annealing of glass and drying 
of textiles. A jet of working fluid leaving a slot or a round 
nozzle is directed towards a targeted heated surface 
where it generates high heat transfer coefficients with 
relatively low pressure drop. The use of slot nozzle pro- 
vides larger stagnation zone with uniform spreading of 
coolant after impingement [1]. In fact, the flow fields and 
heat transfer parameters in slot jets are different from 


Nomenclature 

D, particle diameter (m) 

Cp specific heat capacity of fluid (J/kg K) 
K thermal conductivity (W/m K) 


kp static thermal conductivity (W/m K) 


circular jets [2]. The flow regime can vary from laminar 
to completely turbulent one depending on the application, 
where laminar impinging jets have been used frequently 
for the cooling of electronic chips [3]. Moreover, several 
studies have been conducted on laminar [4] and turbulent 
[5] flow regimes of confined jet impingement. 

Jet Impingement has two major design configurations 
i.e. confined and unconfined. The presence of a confining 
top produces significant effects on fluid flow dynamics 
as well as on heat transfer characteristics of the jet [6]. 
Limited space requirement for compact design makes the 
confined configuration more suitable, while unconfined 
impinging jets are simple in design [7]. The importance 
of confinement in jet impingement heat transfer has been 


U, V, W velocity components 
Xyz spatial coordinates 
Greek symbols 

e shear rate (1/s) 


H nozzle to plate distance (m) 
W jet width (m) 
X length of copper plate (m) 
p pressure (Pa) 
Ap pressure drop (pa) 
q heat flux of the hot wall (W/m’) 
Re Reynolds number 
Pe Peclet number 
Nu Nusselt number 
h heat transfer coefficient (W/m? K) 
T temperature (K) 
Tı lower melting temperature (K) 
T2 upper melting temperature (K) 
Tm melting range (K) 
Lin melting point (K) 
Ay latent heat of fusion of PCM (kJ/kg) 
V velocity 


investigated in [8]. Lin et al. [9] ex-perimentally ob- 
served that heat transfer coefficients increases with the 
jet Reynolds number. Park et al. [10] numerically re- 
ported that the nozzle to target height influence the local 
and average Nusselt numbers for both laminar and tur- 
bulent cases. 

In recent years, researchers mainly focus on different 
types of working fluids used for impinging jet heat 
transfer. Numerous studies has been conducted on air as a 
working fluid for jet impingement [11]. However, im- 
pinging jets with liquid have attracted much more atten- 
tion recently due to its better heat transfer performance 
[12]. Literature suggested that thermophysical properties 
of the working fluid will greatly influence the cooling 
capability of jet impingement. Moreover, the specific 
heat capacity of the traditional fluids like water is not 
sufficient enough to meet the requirements of high heat 
flux removal [13]. To address this issue, an innovative 
technique has been suggested recently to add na- 
no-encapsulated phase change material (NEPCM) par- 
ticles in the base fluid to form a two phase suspension 
[14]. The particle size vary from micro to nano depend- 
ing on the application [15]. 

The PCM has the capability to absorb and release heat 
during melting and solidification simultaneously, while 
encapsulation prevent PCM from leakage [16]. In general, 
NEPCM particles are made up of organic paraffin cores 
with surrounding shells of cross linked polymer [17]. The 
latent heat of NEPCM capsule suspended in the base 
fluid drastically increase the thermal storage capacity of 
the slurry when particles undergoes phase change process 
[18]. Therefore, using these types of advanced fluids are 


Cm mass concentration of NEPCM 


u dynamic viscosity (mPa.s) 

p density (kg/m’°) 

Q viscous dissipation 

a volume concentration of NEPCM 
a thermal diffusivity (m”/s) 
Subscripts 

0 stagnation point 

avg average 

p particle 

b bulk 

w impingement wall 

pcm phase change material 
f fluid 

s solid 
J jet 

eff effective thermophysical properties of fluid 


beneficial for applications such as compact heat ex- 
changer and heat sinks if the cooling system parameters 
are well designed to take maximum advantage of the 
NEPCM latent heat [19]. 

Numerous studies reported the heat transfer enhance- 
ment in forced convection by adding NEPCM particles to 
the coolant. Wu [20] conducted experimental study to 
investigate the effect of NEPCM slurry on heat transfer 
enhancement of jet impingement cooling. The results 
suggested that the volume fraction of the NEPCM in base 
fluid greatly influence the thermal performance. Slurry 
with 28% volume fraction of NEPCM enhance the heat 
transfer coefficient by 50% for jet impingement as com- 
pare to base fluid. 

In the present article, a three dimensional conjugate 
heat transfer model is used to study the cooling and hy- 
drodynamic performance of a confined slot jet impinge- 
ment in laminar regime with NEPCM slurry as a coolant. 
To the best of authors’ knowledge, no published data 
exist in literature on experimental and numerical study of 
confined slot jet impingement using NEPCM slurry. 
Therefore, this attempt is expected to provide a break- 
through in engineering applications such as electronic 
cooling and material processing. 


Geometrical configuration and Model descrip- 
tion 


The schematic diagram of the physical domain is pre- 
sented in Figure 1,which illustrates geometrical configu- 
rations of confined slot jet impingement cooling system. 
The three-dimensional model has nozzle width (W) equal 


to 6.2mm and plate’s length to width ratio (L/W) equal to 
60 while the nozzle to plate distance (H/W) ranging from 
4 to 6. The Cartesian coordinate system is employed in 
the present study, which origin is located at the center of 
impingement surface. A constant heat flux value of 
10,000 W/m? is applied on the bottom surface of copper 
plate. The confinement top is made up of aluminum, 
which is insulated completely so that no heat transfer 
take place across it. Furthermore, uniform inlet velocities 
are considered with laminar and incompressible flow 
regime. Radiation and natural convection effects are 
small enough to be neglected in the model design. 
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Fig. 1 Schematic diagram of confined slot jet impingement 
cooling system 


As present simulation is conducted at low Reynolds 
number, effects caused by viscous dissipation are neg- 
lected [16]. The working fluid is either water or slurry 
with different mass concentrations of NEPCM particles 
in water. Fluid jet temperature is set at 298K for all stu- 
died configurations to ensure that NEPCM particles ex- 
perience phase change process as soon as slurry reaches 
the impingement surface. The particles mass concentra- 
tion in the slurry is below 0.3; therefore, the fluid is as- 
sumed to be Newtonian [21]. Furthermore, it is assumed 
that the base fluid and NEPCM particles flow at the same 
velocity without having any lag between the phases [22]. 
And also, no mass transfer take place between the base 
fluid and NEPCM capsule in order to make it ensure that 
melted PCM inside the capsule will not disperse in the 
base fluid. NEPCM particles are considered to be spher- 
ical in shape with mean diameter of 100nm [18]. Further, 
there is no temperature gradient present inside the par- 
ticle and its melting range is between 293.15K to 
303.15K. The effect of shell material on heat transfer is 
neglected as it is very thin relatively to the core [19]. It is 
also important to mention that, due to homogenous dis- 
tribution of particles across the base fluid, the bulk prop- 
erties of slurry are expected to be constant apart from 
thermal conductivity and specific heat capacity that de- 
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pends on micro-convection and slurry operating temper- 
ature respectively [17]. 


Governing equations 


The continuity, momentum and energy equations are 
considered to be the governing equations of flow and 
heat transfer in our analysis. They are specified by the 
following expressions: 

Continuity equation: 
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Energy equation for the fluid domain: 
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Energy equation for the solid domain: 


a Z2 k, es ilk, a )=0 (6) 
ox\ ` Ox ) Oy Oy) Oz\ a 


The term @ in energy equation predict the rate at 
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which viscous energy is dissipated per unit volume by 
virtue of fluid viscosity. 


Thermophysical properties of NEPCM slurry 


The bulk properties of slurry are function of NEPCM 
concentration and operating temperature. The slurry con- 
sidered in the present analysis is composed of water as a 
base fluid and NEPCM capsules suspended in it. Whereas, 
NEPCM capsules are made up of n-octadecane core as a 
phase change material with melting point (Tm) of 298.15 
K and a surrounding cross linked polymer shell. Table 1 
illustrate thermophysical properties of slurry components. 
The effective properties of slurry are predicted by the 


following co-relations derived in literature: 
Density: 
Pef = (- Cm Water + CmP p (7) 
Where c,, and p, are mass concentration of slurry and 
density of NEPCM particle. 


Table 1 Thermophysical properties of NEPCM slurry com- 
ponents [12,20] 


P Cp k his u 
Octadecane 850 1800 0.34 220.3 - 
Nano Particles 1055 1965 0.22 107.1 - 
Water (292K) 999 4182 0.60 - 0.993 
Water (298K) 996 4143 0.61 - 0.9 
Viscosity [23]: 

-2.5 
Hef = Hwater Gei] (8) 


Where ¢ is the volume fraction or particle loading 
fraction of slurry and Hwater is function of temperature. 

Thermal conductivity: 

Maxwell’s co-relation is used to predict static thermal 
conductivity of NEPCM slurry at rest. It is expressed as 


follows: 
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ky = Bes 


The effects caused due to particle-particle, particle- 
fluid and particle-wall interactions are lumped together to 
compute the effective thermal conductivity of slurry, 
specified by the following co-relation [24]: 


keg = ky (1+ BEPe? | 
B=3, m=1.5, Pe, <0.67 


B=18, m=0.18, 0.67<Pe,<250 0) 
B=3, yo Pe, >250 
11 k 
where Pe, is particle Peclet number which is defined as; 
eD? 
Pe,=—* (11) 
a 
water 


where Qyarey is thermal diffusivity of water and e is 
shear rate, which is a function of spatial coordinates and 
corresponding velocities. Shear rate can be expressed as 
follows: 
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Heat Capacity: 
Cp ef z a — Cm )C p water + CmCp,p (13) 
Where Cp, water and C,, p are specific heat capacity of 


water and NEPCM particle. Alisetti and Roy [25] sug- 
gested that the difference between various profiles for 
calculating the specific heat capacity of NEPCM particle 
is less than 4%. 


C, (kJ/kg.K) 


Fig.2 Sine profile of specific heat capacity model 


Therefore, sine profile is used in the present study to 
predict the specific heat capacity of NEPCM particle as 
illustrated in Figure 2. The melting process starts when the 
temperature of NEPCM particle reaches T, (solidus tem- 
perature) and terminates at T» (liquidus temperature). 
Furthermore, Tm: and Tm are the melting range and melt- 
ing point of NEPCM particle. When the temperature is 
below and above the melting range, specific heat capacity 
of particle is given by Cp, pem. While, for the temperature in 
the melting range; the specific heat capacity is given by: 


h . T-T 
Cpp = Cy pom dafa ~ Cm sina] E] 
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(14) 
More details about co-relation utilization in the present 
study are given in the specified references. 


Numerical procedure and boundary conditions 


A commercial computational fluid dynamics (CFD) 
code FLUENT [26] has been employed in the present 
study as numerical solver. The three-dimensional go- 
verning equations of continuity, momentum and energy 
are discretized by the finite volume method. The convec- 
tion/diffusion terms present in momentum and energy 
equations are discretized with QUICK scheme. Further, 
SIMPLE algorithm is employed in the numerical solver 
to provide a coupling between velocity and pressure. 

The three-dimensional computational grid is displayed 
in Figure 3. A commercial grid generating software 
Gambit is used to divide the computational domain into 
small hexahedral elements in both fluid and solid do- 
mains. A relatively fine grid is adopted in the local do- 
main near the impingement surface, where the solid liq- 
uid interface occurs for conjugate heat transfer boundary 
condition. Coarse grid is used in the rest of domain. 
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Fig.3  Three-dimensional computational grid design 


Numerical solution is terminated when convergence 
has been ensured at a particular iteration. To check if the 
convergence is achieved the summation of absolute val- 
ues of relative errors (i.e. Residuals) of pressure, velocity 
components and temperature falls below 10°, 10° and 
10°, respectively. Furthermore, the Nusselt number and 
temperature values at stagnation point are also checked 
as a second criteria of convergence. 

The assigned boundary conditions are as following: 

Inlet boundary condition: 
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Outlet boundary condition: 
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Bottom wall heat flux boundary condition: 
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Grid independency and model verification 


The grid sensitivity of numerical model was examined 
by considering four grid sizes of 134,872 (coarse), 
369,700 (medium), 876,732 (fine) and 1,798,236 (very 
fine). Grid independency test was conducted for worst 
Reynolds number and water is used as a working fluid. 
The test results predicted that numerical solution become 
grid independed at fine grid and maximum deviation of 
stagnation Nusselt numbers and pressure drops among 
fine and very fine grids is below 3%. Therefor,fine grid 
with 876,732 elements is appropriate enough to capture 
heat and flow characteristics precisely. 

To the best knowledge of the author, no published ex- 
perimental data is present in the literature related to con- 
fined slot jet impingement cooling systems working in 
laminar regime. Therefore, in order to ensure the code 


validity and solution reliability, the numerical simulation 
is conducted with the same geometrical and operating 
parameters as presented in Di Lorenzo [12]. In this vali- 
dation test, water will be taken under consideration as a 
coolant. Furthermore, a constant temperature of 313K is 
maintained at the impingement surface, while all other 
surfaces are insulated in solid domain. The fluid jet tem- 
perature was kept at 292K. Further details regarding the 
numerical models and procedures can be found in the 
original articles [9] and [12]. 

The validation of the present model is illustrated in 
Figure 4 which predicts the Nusselt numbers at stagnation 
point. The maximum deviation found is below 20% which 
show good agreement as Nusselt number calculations are 
extremely sensitive regarding grid size and design. 
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Fig. 4 Verification of present study with past numerical models 
Results and discussion 


Numerical simulation is carried out to evaluate ther- 
mal and hydrodynamic characteristics of confined slot jet 
impingement in laminar regime by using NEPCM slurry 
as coolant. The effects of Reynolds numbers varying 
from 100 to 600 are taken under consideration. The 
ranges of NEPCM particles mass concentration(C,,) and 
nozzle to plate distance (H/W) are C,=0 — 0.28 and 
H/W=4 - 6, respectively.The results are analyzed and 
reported in terms of Nusselt number, heat transfer coeffi- 
cient, pressure drop, and bulk fluid temperature. 

Figure 5 illustrates the effects of mass con-centrations 
of NEPCM particles in slurry on bulk fluid temperature 
at different Reynolds number and constant nozzle to plate 
distance. As seen, using NEPCM slurry as a coolant re- 
duces the bulk temperature of the fluid in comparison 
with water. In addition, increasing the mass concentration 
of NEPCM in slurry reduces the bulk fluid temperature 
considerably. This is due to the fact that adding NEPCM 
capsules to the base fluid increase the effective heat ca- 
pacity of the coolant, which further enhance thermal 
energy storage capability with little rise of bulk fluid tem- 
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perature. An additional observation is that, with increas- 
ing Reynolds number, the fluid bulk temperature de- 
creases and heat transfer enhancement takes place. This 
can be explained with the basic theory that at low Rey- 
nolds number, fluid flows slowly over the impingement 
surface, thus having much more time to absorb heat. 
Therefore, diffusion is the only source for heat transfer 
resulting in higher bulk fluid temperature. Apart from 
this, at higher Reynolds number, associated velocities 
increased and forced convection become dominant factor 
for heat transfer. Hence, in this case more heat will be 
transferred with little rise in the bulk temperature of fluid. 
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Fig. 5 Effect of NEPCM particle concentration on bulk fluid 
temperature 


Figure 6 shows velocity stream lines of impinging 
coolant in three-dimensional domain. The fluid jet shows 
a sudden suppression in velocity just before hitting the 
heated surface. This velocity reduction tends to increase 


the static pressure by converting kinetic energy of par- 
tially stagnated coolant to potential energy. Therefore, 
stagnation zone corresponds to high static pressure. 

The temperature distributions in fluid domain are 
shown in Figure 7 for water and NEPCM slurry with 
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different mass concentrations. As expected, for all the 
cases, thermal boundary layer is relatively thinner at 
stagnation zone as compared to the area near the outlet. 
This is primarily due to the fact that, strong temperature 
gradients exists at stagnation zone as indicated by dense- 
ly packed isotherms. Temperature gradients decreases 
away from the stagnation zone, where thick boundary 
layer is present. 

This means that, higher heat transfer rates will be found 
at stagnation zone. It is also evident from Figure 7 that the 
thermal boundary layer grows slowly across the copper 
plate for NEPCM slurry as compared to water as coolant. 
This is mainly because of latent heat storage of NEPCM 
capsules suspended in water. Another important result 
shown in Figure 7 is that the thermal boundary layer 
thickness reduces as the mass concentration of NEPCM 
increases in base fluid at a fixed Reynolds number and 
nozzle to plate distance. It can be observed that thermal 
layer of water is more thick than the slurry with C,,=0.28. 
Therefore, this reduction in boundary layer thickness will 
enhance heat transfer due to establishment of strong 
temperature gradients across the copper plate in thermal 
boundary layer region. 

Figure 8 depicts the effects of nozzle to plate distance 
(H/W) on stagnation point Nusselt number with jet speed 
variation at a constant NEPCM concentration. This also 
holds true for other mass concentrations of NEPCM par- 
ticles in slurry. It is observed that the Nusselt number at 
stagnation point increases with jet speed at a particular 
nozzle to plate distance. In addition, as H/W increases 
slightly for a fixed jet speed, the stagnation Nusselt 
number decreases considerably. This effect is more pro- 
minent at high jet speeds, for example at 0.1 m/s. The 
common conclusion is that the confined jet impingement 
cooling systems working at low nozzle to plate distance 
is more beneficial for heat transfer enhancement. 


Fig. 6 Slot jet impingement 3D velocity streamlines 
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Fig.7 Effect of NEPCM mass concentration on thermal boundary layer 
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Fig. 8 Effect of Nozzle to plate distance on Nusselt number at 
stagnation point 


In Figure 9, the average heat transfer is plotted against 
the Reynolds number for various slurry concentrations at 
H/W=4. It shows that heat transfer from the impingement 
surface is highly affected by using NEPCM slurry instead 
of water as coolant. Predicted trend shows an enhance- 
ment of heat transfer coefficient with increase of jet 
Reynolds number and NEPCM mass concentration. This 
enhancement is basically a result of high thermal heat 
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Fig. 9 Average heat transfer coefficient as a function of 
NEPCM concentration 


capacity of the slurry. A significant heat transfer im- 
provement is found at NEPCM C,=0.28 as compared 
with water. Another important point to be discussed is 
that the average heat transfer co-efficient increases as 
Reynolds number increases for all considered cases. 

The pressure drop between the nozzle exit and outlet 
of confined slot jet impingement systems for varying 
Reynolds numbers is illustrated in Figure 10 at H/W=4. It 


is obvious that at higher Reynolds number, system exhi- 
bits a significantly higher pressure drops. Moreover, pre- 
ssure drop in system is extremely sensitive with NEPCM 
concentration. This is due to the fact that higher NEPCM 
concentrations in the base fluid results in higher effective 
viscosity of coolant which, in consequence increase 
pressure drops and pumping power of the system. The 
pressure drop recorded for NEPCM C,,=0.28 at Re=600 
is the worst scenario in the present model. Therefore, it is 
highly recommended that the higher pressure drops asso- 
ciated with high concentration of NEPCM particles at 
higher Reynolds number should be considered carefully 
for designing a confined slot jet impingement system 
with NEPCM slurry as coolant. 
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Fig. 10 Effects of Reynold numbers and NEPCM concentra- 
tion on pressure drop 
Conclusion 


The cooling performance and hydrodynamic features 
of NEPCM slurry were compared with that of water, used 
inside the confined slot jet impingement system as coo- 
lants. The effects of jet Reynolds number, jet to plate 
spacing and mass concentration of NEPCM particles on 
heat transfer coefficient, pressure drop and fluid bulk 
temperature were discussed. A good agreement was 
found between the present results and previous numerical 
works. Analysis predicted that adding NEPCM particles 
in the base fluid helps in improving the Nusselt number 
and decrease the bulk temperature of fluid. It was also 
found that the thermal performance of the system highly 
depends on the jet to target spacing. For the Reynolds 
number considered, an increase in mass concentration 
above 0.2 results in very little heat transfer enhancement. 
However, the associated pressure drop will increase dras- 
tically which, in consequence, will reduce the total effi- 
ciency of the system. Therefore, the present findings are 
expected to provide guidelines in designing and optimiz- 
ing the confined slot jet impingement cooling systems 
using NEPCM slurry as a working fluid in the near fu- 
ture. 
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